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ABSTRACT 

Spacecraft for long duration deep space missions will 
need to be designed to survive micrometeoroid 
bombardment of their surfaces some of which may 
actually be punctured. To avoid loss of the entire mission 
the damage due to such punctures must be limited to 
small, localized areas. This is especially true for power 
system radiators, which necessarily feature large surface 
areas to reject heat at relatively low temperature to the 
space environment by thermal radiation. It may be 
intuitively obvious, that if a space radiator is composed of 
a large number of independently operating segments, 
such as heat pipes, a random micrometeoroid puncture 
will result only in the loss of the punctured segment, and 
not the entire radiator. Due to the redundancy achieved 
by independently operating segments, the wall thickness 
and consequently the weight of such segments can be 
drastically reduced. Probability theory is used to estimate 
the magnitude of such weight reductions as the number 
of segments is increased. An analysis of relevant 
parameter values required for minimum mass 
segmented radiators is also included. 

INTRODUCTION 

The literature on space radiator survivability (e.g., 
English and Guentert 1961; Juhasz 1998) confirms 
logical reasoning by indicating that if a space radiator 
were composed of a multiplicity of independently 
operating segments, a random micrometeoroid puncture 
of the radiator would result in the loss of only the 
punctured segment, not the entire radiator. Hence, the 
strategic advantage of segmented radiators for space 
service applications is obvious. Such parallel redundancy 
can be achieved with pumped loop heat rejection 
radiators dividing a single closed circuit carrying cooling 
fluid into a number of independent parallel circuits. But 
the need for separate circulating pumps, valves, and fluid 
ducts makes such a system unwieldy for more than a 
few redundant loops. A much more efficient approach is 
to utilize the inherent parallel redundancy that heat pipes 
offer for space heat rejection applications. Since each 
heat pipe with its extended heat transfer surface fins can 
be considered as a radiator element that is independent 
of all other elements and thus expendable the advantage 


of a heat pipe radiator in terms of modularity and 
minimum complexity is apparent. Since the design of a 
space radiator can accommodate the loss in functionality 
of a certain fraction of elements over the life of a given 
mission, much thinner walls can be specified. As a result 
dramatic savings in weight can be realized when 
compared to a nonsegmented radiator having the same 
overall survival probability. 

Hence there is a clear strategic advantage of segmenting 
large surface area radiators for space applications. The 
purpose of this paper is to use basic probability theory in 
examining the tradeoff between number of segments and 
the segment wall thickness in order to achieve minimum 
radiator mass, subject to certain assumptions regarding 
the total number of segments N s that are required to 
survive till completion of a mission. An analysis of 
requirements for minimum mass systems is also 
included in the appendix to resolve some initial 
discrepancies with results reported in the literature. 

MATHEMATICAL ANALYSIS 

The tradeoff between the number of segments and the 
segment wall thickness, in order to achieve minimum 
radiator weight, can be derived on the basis of a simple 
probability analysis: 

If we assume that for each of N identical radiator 
segments the probability of not being punctured is p, 
then the probability S of having at least N s segments not 
punctured during the spacecraft radiator’s design life is 
given by the cumulative binomial distribution as 


n=N 


S-I 


m 


~n\(N-n)l 


(i -pT n p n 


(i) 


For each segment the probability p of no puncture is 
given by 

p = e~ ha , or \n(\/ p) = hat ; and h = 

( 2 ) 

where h is the average puncture rate per unit time and 
unit area; a is the vulnerable area of each segment; and t 
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( 6 ) 


is the required useful life of the radiator. Note, that for 1/3 

large numbers of segments N, the value of p - N s /N. ^ ~ \^s a ) w \^s at ) 


The puncture rate h is a function of both the depth of 
penetration and the rate of flux of meteoroids. Summers 
and Charters (1959) observed that the wall thickness, T w , 
which a meteoroid will barely penetrate, is proportional to 
the cube root of the meteoroid mass, m. Moreover, 
Whipple (1958) showed that the meteoroid flux (number 
frequency) is inversely proportional to meteoroid size 
(i.e., mass). Assuming that the meteoroid velocity is 
essentially independent of mass, we can express these 
trends as a simple relationship: 

h ~ m _1 ~ T ~ 3 ; or T w = /T 1/3 (3) 


Since the required radiator area at the end of the mission 
and the required mission life are input parameters for the 
analysis and thus are known a priori, the product ( N s a t ) 
can be considered a constant. Segmenting influences 
the weight per unit area (W/N s a) required at the end of 
life by affecting the relative weight term w. The term w, 
that is relative weight per unit segment area, decreases 
continually as the number of segments increases. 

A brief paper on radiator segmenting by English and 
Guentert (1961) reported, without showing assumptions 
or calculation details, that for large values of N the 
minimum weight per unit segment area will occur at 


Using (2) the wall thickness can be written as 


ln(1/p) 


at 


-1/3 


(3.1) 


Multiplying the numerator and denominator of (3.1) by 
(N s )' ,/3 permits rewriting the approximation as 


( /V s^Ln= e ' 1/4 = 0-7788 (7) 

This result differed somewhat with the detailed analysis 
performed by the author (Juhasz 1998) under the 
constraint that a constant number of surviving segments 
N s was required. In this case, the value of NJN for 
minimum weight per unit segment area was found to be 


T w =[N s \n(Vp)Y V3 [N s at] V * (3.2) 

With the objective of minimizing the initial weight of the 
vulnerable radiator surface material while ensuring a 
specified nonpunctured area at the end of radiator design 
life, we can write the following expression: 

W/N s a = pT w N/N s (4) 

Equation (4) can be recast as 

W = (N s a)p(N/N s )r w (4.1) 

where W is the initial weight, p is the density of the 
radiator surface material, and N s a is the required final 
area. In minimizing W/(N S a), for simplicity, only the 
weight of the vulnerable surface area is considered. The 
fin area, manifold ducting, and structure are not taken 
into account. However, since the mass ratio of such 
peripheral components to vulnerable radiator mass can 
be readily determined, one can apply the proper 
correction factor to arrive at an overall radiator mass. 
Also, since material density, p, may be assumed to be a 
constant for a given radiator, it becomes an invariant in 
the analysis. Using (4) and (3.2) a relative weight w can 
then be defined as 

w = (A// N s )[A/ s ln(l/ p)] -1/3 (5) 

Then by substituting (5) and (3.2) into (4.1) the expression 
for total weight shown can be approximated as 


( / V W Ln = e ' 1/3 = 0.7165 (8) 

To resolve the discrepancy the minimum relative mass 
calculations are performed here in detail for the 
alternative constraint of a constant number of initial 
segments N. Two alternate approaches for each of the 
two constraint assumptions are included in the 
Appendix of this report. As shown in section B of the 
appendix, holding the initial number of segments 
N constant, will indeed result in minimum w at N/N = 
e M =0.7788. 

Both assumptions, constant N s or constant N are valid. 
But the case that specifies a required number of 
surviving segments N s is somewhat more conservative, 
since the number of initial segments for a minimum mass 
system would have to be about 8 percent. However, the 
relative mass versus segment ratio curves are fairly flat 
near their minimum values, and both assumptions yield 
results which are approximately within a ±4 percent band 
of 0.75. Hence, for design guidelines, the value of Ns /N 
for minimum weight per unit segment area can be 
assumed as 

( / V /V U~ 0-75 (9) 

Using the equations of this analysis, to solve for S, p, w, 
and and arbitrarily defining the weight and wall 
thickness as unity for a non-segmented radiator with a 
survival probability of 0.999, a computer code was 
developed to tabulate (see Table 1) and plot relative 
weight w and wall thickness T w for values of N (ranging 
from 10 to 1000) as a function of N s IN , as shown in 
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TABLE 1 : COMPARISON OF RESULTS FOR OVERALL SYSTEM RELIABILITY S = 0.999 



fa) Constant N c 

= 717 




(b) Constant N = 

1000 


N s /N 

P 

Relative 

Relative 

N 


N s /N 

P 

Relative 

Relative 

N s 



Weight 

Thickness 





Weight 

Thickness 



w 

T„ 





w 



.30000 

.32928 

.03597 

.01079 

2390.0 


.30000 

.34549 

.04880 

.01464 

300.0 

.40000 

.43591 

.02972 

.01189 

1792.0 


.40000 

.44791 

.03651 

.01460 

400.0 

.50000 

.54038 

.02627 

.01314 

1434.0 


.50000 

.54824 

.02986 

.01493 

500.0 

.60000 

.64273 

.02445 

.01467 

1195.0 


.60000 

.64660 

.02606 

.01564 

600.0 

.68000 

.72296 

.02392 

.01626 

1054.0 


.68000 

.72379 

.02437 

.01657 

680.0 

.69000 

.73264 

.02390 

.01649 

1039.0 


.69000 

.73334 

.02423 

.01672 

690.0 

.70000 

.74257 

.02391 

.01674 

1024.0 


.70000 

.74286 

.0241 1 

.01688 

700.0 

.71000 

.75275 

.02394 

.01700 

1009.0 


.71000 

.75236 

.02401 

.01705 

710.0 

.71500 

.75760 

.02396 

.01713 

1002.0 


.71500 

.75710 

.02396 

.01713 

715.0 

.71700 

.75899 

.02394 

.01717 

1000.0 


.71700 

.75899 

.02394 

.01717 

717.0 

.75000 

.79094 

.02416 

.01812 

956.0 


.75000 

.79008 

.02376 

.01782 

750.0 

.77880 

.81897 

.02455 

.01912 

920.0 


.77880 

.81714 

.02379 

.01853 

779.0 

.80000 

.83867 

.02493 

.01994 

896.0 


.80000 

.83655 

.02392 

.01914 

800.0 

.90000 

.92960 

.02971 

.02674 

796.0 


.90000 

.92621 

.02710 

.02439 

900.0 

1 .0000 

1 .0000 

1 .0000 

1.0000 

717.0 


1.00000 

1.0000 

1 .0000 

1.0000 

1000.0 




Proportion of segments surviving, N s /N 


Figure 1.— Segmented heat pipe radiator characteristics for survival probability S = 0.999. (a) Relative 
weight of vulnerable surface (nonsegmented = 1). (b) Relative wall thickness (nonsegmented = 1). 


figure 1 . Note the dramatic reduction in both weight and wall 
thickness that accrues because of radiator segmenting, 
even for as few as 10 segments, but especially by increasing 
the number of segments N to 1 000. 

To verify the theoretically derived N s IN ratios at which 
minimum relative weight is achieved, table 1 shows 
results for relative weight and thickness for 1,000 
segments (N = 1 ,000). Note that the minimum w value 
(bold number) for constant N s = 717 occurs at 
NJN ~ 0.69, rather than the 0.717 value derived in the 
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appendix. Similarly, for constant N = 1000, the minimum 
occurs at N s IN -0.75 rather than the predicted e m 
value. The reason for this is that the error in the 
assumption inherent in the derivations, namely p - NJN , 
was about 5.5 percent, as can be verified by comparing 
the tabulated values of p and NJN. In table 2 the p 
values are more in agreement for the higher N values (N 
= 3000) and the lower S value (used to reduce the 
computations). Hence the minimum w values (bold 
numbers) in table 2 do occur at NJN values which are in 
close agreement with the theory. As mentioned above, 



TABLE 2: COMPARISON OF RESULTS FOR OVERALL SYSTEM RELIABILITY S = 0.600 
(a) Constant A/,= 2151 * I I (b) Constant A/ = 3000 * 


N s /N 

P 

Relative 

Relative 

N 



Weight 

Thickness 




w 

T. 


.30000 

.30131 

.02431 

.00729 

7170.0 

.40000 

.40164 

.01997 

.00799 

5377.0 

.50000 

.50181 

.01754 

.00877 

4302.0 

.60000 

.60192 

.01619 

.00971 

3585.0 

.68000 

.68197 

.01569 

.01067 

3163.0 

.69000 

.69200 

.01567 

.01081 

3117.0 

.70000 

.70210 

.01565 

.01096 

3072.0 

.71000 

.71203 

.01564 

.01111 

3029.0 

.71500 

.71699 

.01564 

.01118 

3008.0 

.71700 

.71889 

.01564 

.01121 

3000.0 

.75000 

.75184 

.01569 

.01177 

2868.0 

.77880 

.78085 

.01585 

.01234 

2761 .0 

.80000 

.80195 

.01603 

.01282 

2688.0 

.90000 

.90128 

.01831 

.01648 

2390.0 

1 .00000 

.99976 

.12512 

.12512 

2151.0 


N s /N 

P 

Relative 

Relative 

AT 



Weight 

Thickness 



w 

r w 


.30000 

.30198 

.03252 

.00976 

900.0 

.40000 

.40211 

.02427 

.00971 

1200.0 

.50000 

.50215 

.01979 

.00989 

1500.0 

.60000 

.60209 

.01718 

.01031 

1800.0 

.68000 

.68197 

.01597 

.01086 

2040.0 

.69000 

.69195 

.01587 

.01095 

2070.0 

.70000 

.70193 

.01577 

.01104 

2100.0 

.71000 

.71191 

.01569 

.01114 

2130.0 

.71500 

.71689 

.01565 

.01119 

2145.0 

.71700 

.71889 

.01564 

.01121 

2151.0 

.75000 

.75181 

.01546 

.01159 

2250.0 

.77880 

.78072 

.01541 

.01200 

2337.0 

.80000 

.80165 

.01544 

.01235 

2400.0 

.90000 

.90117 

.01697 

.01527 

2700.0 

1.00000 

.99983 

.12512 

.12512 

3000.0 


"Note that due to higher N and lower S, segment reliability p is very close to N/N s as assumed in theory. Hence 
minimum relative weights (bold characters) are for N/N, values as predicted. 


the relative mass versus NJN curves are fairly flat near 
their minimum values, so that the condition shown in (9) 
will be sufficiently accurate for most applications. The 
reader may note, that in table 2 the relative mass value w 
for the nonsegmented radiator ( NJN=1 ) is much lower 
than unity. This is so because the weight of the 
nonsegmented radiator at S = 0.6 is lower than that for a 
nonsegmented one with S = 0.999. The relative weight 
parameter w is referenced to a nonsegmented 
configuration with a survival probability of 0.999. 

Both, table 1 and figure 1 show that for a given number 
of segments N, the relative weight will increase for 
values of N s IN below the value given in equation (9), 
even though the wall thickness decreases. This is 
because as N s decreases, the vulnerable area of each 
segment increases and p, the probability of no puncture 
per segment, decreases. 

As already mentioned in connection with the derivation of 
equation (4), fin area, radiator fluid manifolds, and 
structural weight were not included in this analysis. 
Regarding fin area, the following observation can be 
made: If a certain fraction of the vulnerable area a of 
each segment represents the fin area not rendered 
useless, even if penetrated by a micrometeoroid, then 
the results of the present analysis will be conservative 
(i.e., the relative weights and wall thicknesses will be 
even lower than shown in fig. 1). However, since the 
structure and manifolds will need to be protected by a 
thicker wall or a bumper shield, the benefits provided by 
fins will be offset to some degree, depending on the 
particular design of the spacecraft radiator. 

Nevertheless, the utilization of finned segments that 
increase the radiating area without increasing the 
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vulnerable area will provide an additional margin in the 
design of large segmented space radiators. The use of 
heat pipes, each operating independently with self- 
contained heat transport loops, will further facilitate the 
space radiator design process. 

SEGMENTED RADIATOR DESIGN EXAMPLES 

Having discussed some of the theoretical basis for 
segmenting radiators and for utilizing heat pipes as ideal 
elements or building blocks, let us next discuss a few 
advanced space radiator designs. The conical radiator 
for the 0.1 MW e SP-100 thermoelectric (T/E) space 
power system with a liquid-metal (lithium) cooled reactor 
heat source (Juhasz, 1998) is shown in figure 2. The 
radiator cone is composed of 12 panels, each of which is 
a radiator segment consisting of 226 C-C heat pipes of 
varying condenser length. These heat pipes are 
essentially of the same design as discussed in this 
report. The SP-100 radiator, designed for 7-yr survival 
probability of 0.99, could provide the 2.4 MW, heat 
rejection needs of the power system with only 1 1 of its 
12 panels functioning, each with 75 percent of their 
original number of heat pipes. Thus each of the twelve 
panels could absorb the loss of about 56 heat pipes for a 
total reduction of radiator area equal to that represented 
by over 670 heat pipes, without significant impact on the 
required output power. Note that additional losses in heat 
pipe radiating surface area would result in a rise in T/E 
cold junction temperature beyond the design value at the 
end of life (EOL). As a result the thermal efficiency and 
thus electric power output would be lowered. It should be 
noted, however, that no drastic power system failure 
would occur! Figure 3 shows a proposed heat rejection 
subsystem for a 15-MW a closed-cycle gas turbine power 
system with a nuclear high-temperature gas reactor heat 





Figure 2.— SP-100 advanced heat pipe segmented radiator. 


Inset: Edge view of double C-C face sheet configuration 



Finned heat pipe condensers 



Figure 3.— Flat segmented heat pipe radiator for a nuclear triple loop gas turbine power system. 
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Figure 4.— Heat pipe radiator temperature profiles for a gas turbine power system. (Letters between lines 
designate heat pipe working fluids.) 


source (Juhasz et al., 1993). This segmented radiator 
system is initially composed of 1600 finned heat pipe 
segments of which 1 200 are needed to survive until the 
end of the mission. This radiator was designed to meet 
the heat rejection requirements of three separate gas 
turbine loops. A bifurcated manifold duct, which carries 
the He turbine exhaust stream, communicates thermally 
with 0.05-m diameter, 1.5-m long C-C heat pipe 
evaporators by combined conduction-convection heat 
transfer. The specific weight of such a radiator was 
shown to be less than 1.5 kg/m 2 as compared to 
-15 kg/m 2 for state-of-the-art (SO A) technology. 

The computed temperature profiles for the He gas and 
the heat pipe evaporator wall along the radiator duct are 
shown in figure 4. These profiles were computed under 
the assumption of constant incremental radiator area 
with incremental manifold duct length. To accommodate 
the more than 800-K radiator temperature drop for this 
high-temperature gas turbine power system (turbine inlet 
temperature of 2000 K), the radiator was divided into five 
temperature zones, each of which is defined by the 
operating temperature range associated with a particular 
heat pipe working fluid: sodium, potassium, cesium, 
sulfur-iodine (in lieu of mercury), and water. Note that the 
double C-C face sheet configuration suggested by the 
perspective edge view offers even greater survivability 
and structural rigidity at a somewhat higher weight, 
depending on the thickness of the face sheets. 

CONCLUDING REMARKS 

A derivation of the mathematical relationships required to 
evaluate weight and reliability advantage of radiator 
segmenting was completed, including theoretical 
analyses of conditions for minimum mass systems. A 


computer code was also developed to generate analysis 
results in tabular or graphical form. The drastic reduction 
in mass and wall thickness required to achieve a 
required overall reliability for segmented radiators was 
demonstrated in theory. Several examples of seg- 
mented radiator designs utilizing heat pipe elements 
were described. The advantage of heat pipe technology 
in fabricating large surface area space radiator systems 
with multiple redundancy, which achieve high mission 
reliability at low specific weight, are pointed out in quan- 
titative terms, as compared to non-segmented radiators. 
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APPENDIX 


Determination of Minimum Relative Radiator Mass 


The expression for minimum radiator mass relative to an 
non segmented configuration having a survival 
probability S of 0.999 was shown in (5) as 

w = (A// N s \N s ln(l/ p)]~ 1/3 (A1) 

It was also stated, that for large values of N the value of 
the probability of no puncture of an individual segment, p, 
can be approximated by N s /N, i.e., p = N s /N. 

The analysis of two separate cases is of interest. 


Alternate Method: 

The logarithmic form of (A2) is: 

lnw = Iny- — InA/g --^ln(lny) (A2.0) 

3 3 

Differentiating (A2.0) leads to: 

±*Ul-0-±^ = 0 (A2.1) 

w dy y 3 In y 

Dividing (A2.1) by 1/y it is readily seen that 


A. For a given value of surviving segments N s what value 
of N s /N will minimize w ? 


Letting y = N/N s ~ 1/p 

The relative mass expression w can be written as 

w-y[N s iny]~ 1/3 (A2) 

The minimum of w is found by solving the equation 
represented by dw/dy = 0. 


dw/dy = N~ V 3 W\ 




M/3 *1 


+ 1*[lny]~ 1/3 [- = 0 
(A3) 


Which can be simplified to 
9iv/9y = A/i 1/3 (lny) r4/3 


-l(lnyr 1 + 1 


= 0 


(A4) 


In (A4) the term in brackets yields: 

In y = 1/3 


and, sincel/y = N/N, it follows that N/N = e' 1 3 (A5) 

or (NJN ) = 0.717 (A5.1) 

From the above result one may conclude, that with the 
assumption of a large number of segments, N, for which 
the corollary condition p = N s /N also holds, the relative 
mass w will minimize when the ratio of surviving number 
of segments to initial number of segments is ~ 0.72 if the 
desired number of surviving segments, N s is held 
constant. 


Iny = 1/3 

and N s /N = e- V3 =0.717 (A2.2) 

as obtained in (A5.1) above. 

B. For a given value of initial segments, N, what value of 
N s /N will minimize w ? 

With constant initial number of segments a more 
convenient form of (A1 ) is 

w = {N/N s tN(N s / A/)ln(l/p)] -1/ 3 (B1) 


and with y = N/N s ~ 1/p , wean be written as 

-i-1/3 


w = y 


N 


Iny 


(B2) 


= AT 1 / 3 y 4 / 3 (lnyr 1/3 


(B3) 


Differentiating (B3) and setting the result to zero leads to 


dw/dy = N~ V3 * 


,4/3 


1 -,-4/3 ^ 1 


3 v " yj 3 
which can be simplified to 


+ 4y 1/3 * (Iny) 1/3 1 = 0 


(B4) 


AT 1 / 3 y 1/3 


[i (lnyr 


4/3 f — Iny * (lny)“ 4/3 


(B5) 


and 


/V" 1/3 y 1/3 (lny) _4/3 


1 4 i 

— + — Iny 
3 3 


= 0 


(B6) 
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In (B 6 ) the term in brackets yields 



or y = e 


1/4 


(B7) 


1 A/c No MA 

and since — = — , it follows that — ^ = e - 0.7788 

y N N 


(B 8 ) 


From the result shown in (B 8 ) one may conclude, that 
with the assumption of a large number of segments, N, 
with the corollary condition p ~ N s A N , the relative mass 
w will minimize when the ratio of surviving number of 
segments to initial number of segments is -0.78 if the 
desired number of initial segments N is held constant. 


Alternate Method: Writing the logarithmic form of (B 3 ), 

we have lnw = -— \nN + — Iny- — ln(lny), which can be 
3 3 3 

differentiated to yield 


1 dw 4 1 1 1 1 1 ( A ^ ) 

w dy 3 y 3 Iny y 3 y[ In yj 


(B9) 


which results in Iny = — ; y = e 1/4 and 


No 1 \l A 

— 7 - =- = e~ v 4 = 0.7788 
N y 


(BIO) 


as obtained in (B 8 ) above. 
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